Suppression of Density Fluctuations in a Quantum Degenerate Fermi Gas 
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We study density profiles of an ideal Fermi gas and observe Pauli suppression of density fluctua- 
tions (atom shot noise) for cold clouds deep in the quantum degenerate regime. Strong suppression 
is observed for probe volumes containing more than 10,000 atoms. Measuring the level of sup- 
pression provides sensitive thermometry at low temperatures. After this method of sensitive noise 
measurements has been validated with an ideal Fermi gas, it can now be applied to characterize 
phase transitions in strongly correlated many-body systems. 
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Systems of fermions obey the Pauli exclusion principle. 
Processes that would require two fermions to occupy the 
same quantum state are suppressed. In recent years, sev- 
eral classic experiments have directly observed manifesta- 
tions of Pauli suppression in Fermi gases. Antibunching 
and the suppression of noise correlations are a direct con- 
sequence of the forbidden double occupancy of a quantum 
state. Such experiments were carried out for electrons [lj- 
13], neutral atoms [J, |5|, and neutrons [6(. In principle, 
such experiments can be done with fermions at any tem- 
perature, but in practice low temperatures increase the 
signal. A second class of (two-body) Pauli suppression 
effects, the suppression of collisions, requires a tempera- 
ture low enough such that the de Broglie wavelength of 
the fermions becomes larger than the range of the inter- 
atomic potential and p-wave collisions freeze out. Exper- 
iments observed the suppression of elastic collisions \AM 
and of clock shifts in radio frequency spectroscopy [9lllCj. 

Here we report on the observation of Pauli suppression 
of density fluctuations, a many-body phenomenon which 
occurs only at even lower temperatures in the quantum 
degenerate regime, where the Fermi gas is cooled below 
the Fermi temperature and the low lying quantum states 
are occupied with probabilities close to one. In contrast, 
an ideal Bose gas close to quantum degeneracy shows 
enhanced density fluctuations 
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The development of a technique to sensitively measure 
density fluctuations was motivated by the connection be- 
tween density fluctuations and compressibility through 
the fluctuation dissipation theorem. In this paper, we 
validate our technique for determining the compressibil- 
ity by applying it to the ideal Fermi gas. In future work, 
it could be extended to interesting many-body phases in 
optical lattices which are distinguished by their incom- 
pressibility |12 |. These include the band insulator, Mott 
insulator, and also the antiferromagnet for which spin 
fluctuations, i.e. fluctuations of the difference in density 
between the two spin states, are suppressed. 

Until now, sub-Poissonian number fluctuations of ul- 
tracold atoms have been observed only for small clouds of 
bosons with typically a few hundred atoms }13l - ll6l ] and 
directly [171, |l8( or indirectly [19( for the bosonic Mott in- 




FIG. 1: Phase space diagram of ballistic expansion of a har- 
monically trapped Fermi gas. Ballistic expansion conserves 
phase space density and shears the initially occupied spheri- 
cal area into an ellipse. In the center of the cloud, the local 
Fermi momentum and the sharpness of the Fermi distribution 
are scaled by the same factor, keeping the ratio of local tem- 
perature to Fermi energy constant. The same is true for all 
points in the expanded cloud relative to their corresponding 
unsealed in-trap points. 



sulator in optical lattices. For fermions in optical lattices, 
the crossover to an incompressible Mott insulator phase 
was inferred from the fraction of double occupations [201 ] 



or the cloud size [2l|. Here we report the observation of 



density fluctuations in a large cloud of fermions, showing 
sub-Poissonian statistics for atom numbers in excess of 
10,000 per probe volume. 

The basic concept of the experiment is to repeat- 
edly produce cold gas clouds and then count the num- 
ber of atoms in a small probe volume within the ex- 
tended cloud. Many iterations allow us to determine 
the average atom number N in the probe volume and 
its variance (AN) 2 . For independent particles, one ex- 
pects Poisson statistics, i.e. (AN) 2 / (N) — 1. This is di- 
rectly obtained from the fluctuation dissipation theorem 
(AN) 2 / (N) — nksTKr, where n is the density of the gas, 
and kt the isothermal compressibility. For an ideal clas- 
sical gas kt = l/(nfcflT), and one retrieves Poissonian 
statistics. For an ideal Fermi gas close to zero temper- 
ature with Fermi energy Ep, Kp = 3/(2nEp), and the 
variance (AN) 2 is suppressed below Poissonian fluctua- 
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FIG. 2: (Color online) Comparison of density images to vari- 
ance images. For Poissonian fluctuations, the two images at a 
given temperature should be identical. The variance images 
were obtained by determining the local density fluctuations 
from a set of 85 images taken under identical conditions, (a) 
Two dimensional image of optical density of an ideal Fermi 
gas after 7 ms of ballistic expansion. The noise data were 
taken by limiting the field of view to the dashed region of 
interest, allowing for faster image acquisition, (b) For the 
heated sample, variance and density pictures are almost iden- 
tical, implying only modest deviation from Poissonian statis- 
tics, (c) Fermi suppression of density fluctuations deep in the 
quantum degenerate regime manifests itself through the dif- 
ference between density and variance picture. Especially in 
the center of the cloud, there is a large suppression of density 
fluctuations. The variance images were smoothed over 6x6 
bins. The width of images (b) and (c) is 2 mm. 



tions by the Pauli suppression factor 3kBT/(2Ep). All 
number fluctuations are thermal, as indicated by the pro- 
portionality of (AN) 2 to the temperature in the fluctua- 
tion dissipation theorem. Only for the ideal classical gas, 
where the compressibility diverges as 1/T, one obtains 
Poissonian fluctuations even at zero temperature. 

The counting of atoms in a probe volume can be done 
while the atoms are trapped, or after ballistic expansion. 
Ballistic expansion maintains the phase space density and 
therefore the occupation statistics. Consequently, den- 
sity fluctuations are exactly rescaled in space by the bal- 
listic expansion factors as illustrated in FigQ] [24J, [25]. 
Note that this rescaling is a unique property of the har- 
monic oscillator potential, so future work on density fluc- 
tuations in optical lattices must employ in-trap imaging. 
For the present work, we chose ballistic expansion. This 
choice increases the number of fully resolved bins due to 
optical resolution and dcpth-of-ficld, it allows adjusting 
the optimum optical density by choosing an appropri- 
ate expansion time, and it avoids image artifacts at high 
magnification. 

We first present our main results, and then discuss im- 
portant aspects of sample preparation, calibration of ab- 
sorption cross section, data analysis and corrections for 
photon shot noise. Fig. [2^l shows an absorption image of 
an expanding cloud of fermionic atoms. The probe vol- 
ume, in which the number of atoms is counted, is chosen 
to be 26 fxm in the transverse directions, and extends 
through the entire cloud in the direction of the line of 
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FIG. 3: Comparison of observed variances (black dots) with a 
theoretical model (black line) and the observed atom number 
(gray), at three different temperatures (a, b, and c), showing 
50, 40, and 15% suppression. Noise thermometry is imple- 
mented by fitting the observed fluctuations, resulting in tem- 
peratures T/T F of 0.23±.01, 0.33±.02, and 0.60±.02. This is 
in good agreement with temperatures 0.21±.01, 0.31±.01, and 
0.6±.l obtained by fitting the shape of the expanded cloud 
[29] . The quoted uncertainties correspond to one standard 
deviation and are purely statistical. 



sight. The large transverse size completely avoids av- 
eraging of fluctuations due to finite optical resolution. 
From 85 such images, after careful normalization [26], 
the variance in the measured atom number is determined 
as a function of position. After subtracting the photon 
shot noise contribution, a 2D image of the atom number 
variance (AJV) 2 is obtained. For a Poissonian sample 
(with no suppression of fluctuations), this image would 
be identical to an absorption image showing the number 
of atoms per probe volume. This is close to the situa- 
tion for the hottest cloud (the temperature was limited 
by the trap depth), whereas the colder clouds show a 
distinct suppression of the atom number variance, espe- 
cially in the center of the cloud where the local T/Tp is 
smallest . 

In Fig. [31 profiles of the variance are compared to 
theoretical predictions 
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wavevector q are proportional to the structure factor 
S(q,T). Since our probe volume (transverse size 26 jjxci) 
is much larger than the inverse Fermi wavevector of the 
expanded cloud (1/<?f = 1.1/mi), S(q = 0,T) has been 
integrated along the line of sight for comparison with the 
experimental profiles. Within the local density approxi- 
mation, S(q = 0, T) at a given position in the trap is the 
binomial variance n/j(l— rik) integrated over all momenta, 
where the occupation probability Hk(k, /J-,T) is obtained 
from the Fermi-Dirac distribution with a local chemical 
potential /x determined by the shape of the trap. Fig. [4] 
shows the dependence of the atom number variance on 
atom number for the hot and cold clouds. A statistical 
analysis of the data used in the figure is in [26] . 

The experiments were carried out with typically 2.5 
million 6 Li atoms per spin state confined in a round 
crossed dipole trap with radial and axial trap frequencies 
uj r = 2ir x 160 s _1 and ui z — 2tt x 230 s _1 corresponding 



to an in-trap Fermi energy of Ep — ks X 2.15 /iK. The 
sample was prepared by laser cooling followed by sympa- 
thetic cooling with 23 Na in a magnetic trap. 6 Li atoms 
in the highest hyperfine state were transferred into the 
optical trap, and an equal mixture of atoms in the low- 
est two hyperfine states was produced. The sample was 
then evaporatively cooled by ramping down the optical 
trapping potential at a magnetic bias field B = 320±5 
G where a scattering length of -300 Bohr radii ensured 
efficient evaporation. Finally, the magnetic field was in- 
creased to B = 520±5 G, near the zero crossing of the 
scattering length, realizing a non-interacting Fermi gas. 
Absorption images were taken after 7 ms of ballistic ex- 
pansion. 

We were careful to prepare samples at different tem- 
peratures with similar cloud sizes and central optical den- 
sities to make sure that they were imaged with the same 
effective cross section and resolution. Hotter clouds were 
prepared by heating the colder cloud with parametric 
modulation of the trapping potential. For the hottest 
cloud this was done near 520 G to avoid excessive evap- 
oration losses. 

Atomic shot noise dominates over photon shot noise 
only if each atom absorbs several photons. As a result, 
the absorption images were taken using the cycling tran- 
sition to the lowest lying branch of the 2 P^/2 manifold. 
However, the number of absorbed photons that could 
be tolerated was severely limited by the acceleration of 
the atoms by the photon recoil, which Doppler shifts the 
atoms out of resonance. Consequently, the effective ab- 
sorption cross section depends on the probe laser inten- 
sity and duration. To limit the need for nonlinear nor- 
malization procedures, we chose a probe laser intensity 
corresponding to an average of only 6 absorbed photons 
per atom during 4 /is of exposure time. At this intensity, 
about 12% of the 6 Li saturation intensity, the measured 
optical density was found to be reduced by 20% from its 
low-intensity value |26j |. For each bin, the atom number 
variance is obtained by subtracting the known photon 
shot noise from the variance in the optical density [26J. 

The absorption cross section is a crucial quantity in 
the conversion rate between the optical density and the 
number of detected atoms. For the cycling transition, 
the resonant absorption cross section is 2.14 x 10~ 13 m . 
Applying the measured 20% reduction mentioned above 
leads to a value of 1.71 x 10~ 13 m 2 . This is an upper limit 
to the cross section due to imperfections in polarization 
and residual line broadening. An independent estimate 
of the effective cross section of 1.48 x 10 ~ 13 m 2 was ob- 
tained by comparing the integrated optical density to the 
number of fermions necessary to fill up the trap to the 
chemical potential. The value of the chemical potential 
was obtained from fits to the ballistic expansion pictures 
that allowed independent determination of the absolute 
temperature and the fugacity of the gas. We could not 
precisely assess the accuracy of this value of the cross 
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FIG. 4: Atom number variance vs. average atom number. 
For each spatial position, the average atom number per bin 
and its variance were determined using 85 images. The filled 
and open circles in the figure are averages of different spatial 
bin positions with similar average atom number. For a hot 
cloud at T/Tf—0.6 (filled circles), the atom number variance 
is equal to the average atom number (dotted line, full Pois- 
sonian noise) in the spatial wings where the atom number is 
low. The deviation from the linear slope for a cold cloud at 
T/Tf—0.21 (open circles) is due to Pauli suppression of den- 
sity fluctuations. There is also some suppression at the center 
of the hot cloud, where the atom number is high. The solid 
and dashed lines are quadratic fits for the hot and cold clouds 
to guide the eye. 



section, since we did not fully characterize the effect of a 
weak residual magnetic field curvature on trapping and 
on the ballistic expansion. The most accurate value for 
the effective cross section was determined from the ob- 
served atom shot noise itself. The atom shot noise in the 
wings of the hottest cloud is Poissonian, and this condi- 
tion determines the absorption cross section. Requiring 
that the slope of variance of the atom number (AJV) 2 vs. 
atom number N is unity (see Fig. 2]) results in a value of 
(1.50 ± 0.12) x 10~ 13 m 2 for the effective cross section in 
good agreement with the two above estimates. 

The spatial volume for the atom counting needs to be 
larger than the optical resolution. For smaller bin sizes 
(i.e. small counting volumes), the noise is reduced since 
the finite spatial resolution and depth of field blur the 
absorption signal. In our setup, the smallest bin size 
without blurring was determined by the depth of field, 
since the size of the expanded cloud was larger than the 
depth of field associated with the diffraction limit of our 
optical system. We determined the effective optical res- 
olution by binning the absorption data over more and 
more pixels of the CCD camera, and determining the 
normalized central variance (A N) 2 /N vs. bin size [26j . 
The normalized variance increased and saturated for bin 
sizes larger than 26 /xm (in the object plane), and this 
bin size was used in the data analysis. We observe the 
same suppression ratios for bin sizes as large as 40 /xm, 
corresponding to more than 10,000 atoms per bin. 



For a cold fermion cloud, the zero temperature struc- 
ture factor S(q) becomes unity for q > 2qp. This reflects 
the fact that momentum transfer above 2qp to any parti- 
cle will not be Pauli suppressed by occupation of the final 
state. In principle, this can be observed by using bin sizes 
smaller than the Fermi wavelength, or by Fourier trans- 
forming the spatial noise images. For large values of q, 
Pauli suppression of density fluctuations should disap- 
pear, and the noise should be Poissonian. However, our 
imaging system loses its contrast before q « 2qp [26j. 

Observation of density fluctuations, through the 
fluctuation-dissipation theorem, determines the product 
of temperature and compressibility. It provides an ab- 
solute thermometer, as demonstrated in Fig. [3] if the 
compressibility is known or is experimentally determined 
from the shape of the density profile of the trapped 
cloud [171, |33j . Because variance is proportional to tem- 
perature for T <C Tp, noise thermometry maintains its 
sensitivity at very low temperature, in contrast to the 
standard technique of fitting spatial profiles. 

Density fluctuations lead to Rayleigh scattering of 
light. The differential cross section for scattering light of 
wavevector k by an angle 6 is proportional to the struc- 
ture factor S(q), where q = 2ksin(9/2) [28[. In this work, 
we have directly observed the Pauli suppression of density 
fluctuations and therefore S(q) < 1, which implies sup- 
pression of light scattering at small angles (corresponding 
to values of q inversely proportional to our bin size) . How 
are the absorption images affected by the suppression of 
light scattering? Since the photon recoil was larger that 
the Fermi momentum of the expanded cloud, large-angle 
light scattering is not suppressed. For the parameters 
of our experiment, we estimate that the absorption cross 
section at the center of a T = Fermi cloud is reduced by 
only 0.3% due to Pauli blocking [3l(. Although we have 
not directly observed the Pauli suppression of light scat- 
tering, which has been discussed for over 20 years [30h32i |. 
by observing reduced density fluctuations we have seen 
the underlying mechanism for suppression of light scat- 
tering. 

In conclusion, we have established a sensitive technique 
for determining atomic shot noise and observed the sup- 
pression of density fluctuations in a quantum degenerate 
ideal Fermi gas. This technique is promising for ther- 
mometry of strongly correlated many-body systems and 
for observing phase-transitions or cross-overs to incom- 
pressible quantum phases. 

We acknowledge Joseph Thywissen and Markus 
Greiner for useful discussions. This work was sup- 
ported by NSF and the Office of Naval Research, AFOSR 
(through the Multidisciplinary University Research Ini- 
tiative program), and under Army Research Office grant 
no. W911NF-07-1-0493 with funds from the Defense Ad- 
vanced Research Projects Agency Optical Lattice Emu- 
lator program. 



[2] M. Henny et al., Science 284, 296 (1999). 

[3] H. Kiesel, A. Renz, and F. Hasselbach, Nature 418, 392 

(2002). 
[4] T. Rom et al., Nature 444, 733 (2006). 
[5] T. Jeltes et al., Nature 445, 402 (2007). 
[6] M. Iannuzzi, A. Orecchini, F. Sacchetti, P. Facchi, and 

S. Pascazio, Phys. Rev. Lett. 96, 080402 (2006). 
[7] B. DeMarco, J. L. Bohn, J. P. Burke, M. Holland, and 

D. S. Jin, Phys. Rev. Lett. 82, 4208 (1999). 
[8] B. DeMarco, S. B. Papp, and D. S. Jin, Phys. Rev. Lett. 

86, 5409 (2001). 
[9] M. W. Zwierlein, Z. Hadzibabic, S. Gupta, and W. Ket- 
terle, Phys. Rev. Lett. 91, 250404 (2003). 

[10] S. Gupta et al, Science 300, 1723 (2003). 

[11] J. Esteve et al, Phys. Rev. Lett. 96, 130403 (2006). 

[12] Q. Zhou, Y. K ato, N. Kawashima, and N. Trivedi, 
larXiv:0901.0606l . 

[13] G S. Chuu et al, Phys. Rev. Lett. 95, 260403 (2005). 

[14] J. Esteve, G Gross, A. Weller, S. Giovanazzi, and M. K. 
Oberthaler, Nature 455, 1216 (2008). 

[15] S. Whitlock, C. F. Ockeloen, and R. J. G Spreeuw, Phys. 
Rev. Lett. 104, 120402 (2010). 

[16] A. Itah et al., Phys. Rev. Lett. 104, 113001 (2010). 

[17] N. Gemelke, X. Zhang, C. L. Hung, and G Chin, Nature 
460, 995 (2009). 

[18] M. Greiner, unpublished . 

[19] M. Greiner, O. Mandel, T. W. Hansen, and I. Bloch, 
Nature 419, 51 (2002). 

[20] R. Jordens, N. Strohmaier, K. Giinter, H. Moritz, and 
T. Esslinger, Nature 415, 204-207 (2008). 

[21] U. Schneider et al, Science 322, 1520 (2008). 

[22] M. Greiner, C. A. Regal, J. T. Stewart, and D. S. Jin, 
Phys. Rev. Lett. 94, 110401 (2005). 

[23] S. Foiling et al, Nature 434, 481 (2005). 

[24] S. Gupta, Z. Hadzibabic, J. R. Anglin, and W. Ketterle, 
Phys. Rev. Lett. 92, 100401 (2004). 

[25] G M. Bruun and C. W. Clark, Phys. Rev. A 61, 
061601(R) (2000). 

[26] See EPAPS Document No. [number will be inserted by 
publisher] for additional details. For more information on 
EPAPS, see http://www.aip.org/pubservs/epaps.html 

[27] Y. Castin, in Proceedings of the International School of 
Physics Enrico Fermi, Course CLXIV, edited byM. In- 
guscio, W. Ketterle, and C. Salomon (IOS, Amsterdam), 
2008. 

[28] D. Pines and P. Nozieres, The Theory of Quantum Liq- 
uids, (Addison- Wesley, MA, 1988), Vol. 1. 

[29] W. Ketterle and M. W. Zwierlein, in Proceedings of 
the International School of Physics Enrico Fermi, Course 
CLXIV, edited by M. Inguscio, W. Ketterle, and C. Sa- 
lomon (IOS, Amsterdam), 2008. 

[30] K. Helmerson, M. Xiao, and D. E. Pritchard, in Inter- 
national Quantum Electronics Conference 1990, Book of 
Abstracts, 1990, QTHH4. 

[31] A. Gorlitz, A. P. Chikkatur, and W. Ketterle, Physical 
Review A 63, 041601(R) (2001). 

[32] B. Shuve and J. H. Thywissen, J. Phys. B 43, 015301 
(2010). (and references therein) 

[33] Q. Zhou and T. L. Ho, larXiv:0908. 30151 . 



[1] W. D. Oliver, J. Kim, R. C. Liu, and Y. Yamamoto, 
Science 284, 299 (1999). 



